arXiv: 1507.08915v2 [gr-qc] 23 Jan 2016 


On the stability of dyons and dyonic black holes in 
Einstein-Yang-Mills theory 


Brien C. Nolan^ and Elizabeth Winstanley^ 

^ School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, 
Ireland 

^ Consortium for Fundamental Physics, School of Mathematics and Statistics, 

The University of Sheffield, Hicks Building, Hounsfield Road, Sheffield. S3 7RH 
United Kingdom 

E-mail: brien.nolan@dcu.ie, E.Winstanley@sheffield.ac.uk 
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metric and non-Abelian gauge field. 
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1. Introduction 

Since the discovery of non-trivial soliton pQ and black hole [2] solutions of the four¬ 
dimensional su(2) Einstein-Yang-Mills (EYM) equations in asymptotically flat space- 
time, the EYM system has been studied extensively (see [3] for a review). For 5u(2) 
gauge group, the gauge held of non-trivial solutions in four-dimensional asymptotically 
hat space-time is purely magnetic [1] and, furthermore, these solutions are unstable 
under linear, spherically symmetric, perturbations [5]. 

The properties of EYM solutions in asymptotically anti-de Sitter (adS) space-time 
are very diherent from those in asymptotically hat space-time. The hrst diherence is 
the existence of four-dimensional, spherically symmetric, purely magnetic soliton and 
black hole [7] solutions of su(2) EYM in adS which are stable under linear, spherically 
symmetric, perturbations. Subsequently it was shown that there exist both soliton and 
black hole solutions which are stable under general linear perturbations of the metric 
and gauge held [8]. If the gauge group is enlarged to su(Y), purely magnetic, spherically 
symmetric, soliton and black hole solutions with iV — 1 gauge held degrees of freedom 
exist [9] . It can be proven that at least some of these are stable under linear, spherically 
symmetric perturbations nm. 
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The second surprising feature of solutions of su(2) EYM in adS is the existence 
of non-trivial, spherically symmetric, dyonic solitons and black holes mm- For these 
solutions the gauge held has a non-trivial electric part as well as a magnetic part. 
Properties of these spherically symmetric dyonic solitons and black holes were explored 
numerically in [6]. The existence of non-trivial dyonic solutions in a neighbourhood of 
the trivial (embedded Schwarzschild-adS) solution was proven in |TT]. Although these 
dyonic solutions were discovered numerically over hfteen years ago, their stability has 
remained an open question which we address in this paper. 

We consider static, spherically symmetric, dyonic soliton and black hole solutions of 
su(2) EYM in adS. In section [2] we introduce the action and held equations, and briehy 
review some of the properties of the static equilibrium solutions mm- Next, in section 
m we derive the equations governing time-dependent, linear, spherically symmetric 
perturbations of the static equilibrium solutions. The analysis results in a pair of 
coupled Schrodinger-like equations for two of the perturbations. The third independent 
perturbation is governed by a constraint equation which does not involve any derivatives 
with respect to time. Section 0] contains our proof of the existence of non-trivial dyonic 
solitons and black holes, in a neighbourhood of the embedded trivial solution, which are 
stable under the linear perturbations. Finally we present our conclusions in section |5l 

2. Dyons and dyonic black holes in su(2) Einstein-Yang-Mills theory 

In this section we introduce the action and held equations for su(2) Einstein-Yang- 
Mills theory with a negative cosmological constant. We also briehy review some of the 
properties of the static, spherically symmetric, dyon and dyonic black hole solutions of 
this theory, which were found numerically in [6] and whose existence was proven in m- 

2.1. Ansatz and field equations 

We begin with the action for Einstein-Yang-Mills theory in four-dimensional space-time 
with a cosmological constant A: 

^EYM = lj [R-2A- Tr , (2.1) 

where R is the Ricci scalar, g is the metric determinant and iRj, is the Yang-Mills gauge 
held. Tr denotes a Lie algebra trace. Here and throughout this paper, the space-time has 
signature (—, -f, -f, -|-), we use units in which AnG = 1 = c and we have hxed the gauge 
coupling constant to be equal to unity. We consider a negative cosmological constant 
A < 0 and the gauge Lie algebra is su(2). 

Varying the action (12.11) we obtain the held equations 

‘2Tti/ Rtu T Ag.j-jy^ 


(2.2a) 

{2.2b) 
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where the Yang-Mills held strength tensor takes the form 

Ftu = drAy — dyA^ + [At, Af \, (2-3) 

with Ar the Yang-Mills gange potential, and [Ar, Af\ denoting the Lie algebra 
commntator. The stress energy tensor is 

Ttu = Tr FrxF/ - Tr . (2.4) 

In this paper we are interested in the stability of static, spherically symmetric, dyon 
and dyonic black hole solntions of the held eqnations fl2.2af[Y2^ . In the next section we 
shall consider time-dependent, linear, spherically symmetric pertnrbations of the static 
eqnilibrinm solntions, so we consider a time-dependent, spherically symmetric metric as 
follows 


=—fx{t,r)S{t,rY dt‘^ + fi{t,r) ^ dr'^ + r‘^ [d9‘^ + sin^ 9 d(j)‘^) , (2.5) 


where the metric fnnctions yi{t, r) and S{t, r) depend on time t and the radial co-ordinate 
r. We may write the metric fnnction /i(t, r) in the alternative form 


fx{t,r) = 1 


2 m(t, r) 


+ 


r & ’ 

where the adS radius of curvature i is given by 

^-4 


( 2 . 6 ) 


(2.7) 


The time-dependent, spherically symmetric su(2) Yang-Mills gauge potential At 
can be written as follows, after an appropriate choice of gauge | 12 j : 

A = Adt + Bdr + ^{C -C^) d9-^- [{C + C^) sm9 + D cos9] d^,{2.8) 
where A, B, C and D are 2 x 2 matrices, given by 


A^i( 0 A 

2 0 -a{t,r) ) ’ 


0 oj{t,r)F'^Fr) \ 

0 0 J 


? / I3{t,r) 0 A 

2 0 -/3(t,r) ) 



(2.9) 


Here, a{t,r), (3{t,r), '^{t,r) and oj{t,r) are real functions of time t and the radial co¬ 
ordinate r. The matrix is the Hermitian conjugate of the matrix C. 


2.2. Static, spherically symmetric, dyons and dyonic black holes 

For static equilibrium solutions of the held equations, the metric functions m = mo(r) 
and S = 5'o(r) now depend only on the radial co-ordinate r. By a choice of gauge [T 2 ] . 
the gauge held function ft can be set to zero, and then one of the Yang-Mills equations 
reduces to 7 = 0. The remaining gauge held functions, a = ao{r) and u = Uo{r), are 
also functions of r only. 
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The field equations fl 2 . 2 d \2.2b^ then reduce to the following static field equations: 


2^/2 


r a, 


2, ,2 


rUn = 


0 




2\2 


2^2 


ho 


So 




r 


0 — hoQ^o + 


2/io _ ho>S'o 


(Uq 


200 ^^ 


0 — fiou'f + ^/i'o + '^0 + (l “ '^o) + 


tto^o 


ho5, 


2 ’ 


(2.10a) 

(2.106) 

( 2 . 10 c) 

(2.10d) 


where a prime ' denotes differentiation with respect to r. The static field equations 
fl2.10afl2.10dp possess the following symmetries. Firstly they are invariant under the 
transformation ao —)■ —cto; secondly the transformation ojo —)■ —oJo also leaves them 
unchanged; and finally they are preserved by the scaling symmetry: 


So —^ •^5*0, Oio —y -^cTo, (2-11) 

for any constant A. The scaling symmetry fl 2 . 1 ip arises due to the invariance of the 
static metric and gauge potential under rescalings of the time co-ordinate t —)■ X~^t. 
When the metric and gauge potential are time-dependent, the gauge freedom remaining 
in rescaling the time co-ordinate is discussed in section 13.21 

The static field equations (l2.10af[T10dl) have three singular points of interest. These 
are located at the origin r = 0 (relevant only for soliton solutions), at the black hole 
horizon (corresponding to zeros r = Vh of the metric function /i, if there are any), and 
as r —)■ oo. As pointed out in m, while zeros of the metric function S yield a fourth 
possible singular point, these are not of relevance to the classes of solutions considered 
here. Suitable boundary conditions therefore have to be imposed on the field variables 
at the singular points r = 0,r = and r —)■ oo. Near the origin, the field variables take 
the form mm 

mo{r) = + 2uji^ + O(r^), 

So{i^) = 5*1 -f -f 4S'ia;2^ -f O(r^), 

aoir) = air + y + Suj + 2 u 2 - + O(r^), 

oooir) = l + uj 2 r^+ 0{r^), ( 2 . 12 ) 

and the solutions are parameterized by the constants Si, ai, U 2 and the adS radius 
of curvature i. In a neighbourhood of the black hole event horizon, the corresponding 
expansion of the field variables is mm 

3 

rno{r) = ^ + ^ + m'^{r - Th) + 0{r - rtf, 

S'o(r) = Sh + SUr- r,) + 0{r - 
ao(r) = a(^ (r - Vh) + 0{r - ruf, 
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= ^h + ^'hi.'T - rh) + 0{r - Thf, (2.13) 

where m'f^, S'f^ and u'j^ are given by the field eqnations fl2.10al - [Z10dj) in terms of Sh, 
and Uh- 




2^2 


, _ , (1 - 

^ “r 


2Sl 

^ SkPk^i'irk)'^ 


2rl 


12 


+ 


OJu = 








2ShUj'f 

rh 

1 )> 


where 


/ / X 1 2m'f, 3rh 

/i' r, =-h + ^>o, 

rh rh 


(2.14) 


(2.15) 


so that the constants Sh, Oh and ojh, together with the adS radius of curvature i, 
parameterize the solutions. At infinity the field variables have the following behaviour 

nm- 


mo{r) = M ~ 


9 2 /72 ^1 

y + “ 


' uj: 




+ Oir-^), 


S'o(r) = 1 - — + cf) + 0{r ®), 


2r4 
di 


aoir) = + — + 0{r ^), 

r 

Uci(f) =iPco + — + 0(r~^), 


(2.16) 


where M, ctoo, cjoo, Ci and di are arbitrary constants. The fact that 5'o(r) ^ 1 as 
r —>■ cx) fixes the parameters Si and Sh in the expansions of So near the origin fl2.12p 
and event horizon fl2.13p . respectively. In practice, however, we can regard Si and Sh 
as free parameters, since, if Soo 7^ 1, a scaling transformation fl2.1ip with A = S^^ can 
always be applied. 

The field equations fl2.10af[T10dp possess a trivial solution given by 

ao(r) = 0, Uo(r) = ±1, mo(r) = M, S'o(r) = 1. (2.17) 

For M > 0 this is the Schwarzschild-adS black hole; for M = 0 this is pure adS 
space-time. There are also embedded (electrically and magnetically charged) Abelian 
Reissner-Nordstrom-adS solutions of the static field equations [m, but we shall not 
consider these further in this paper. Purely magnetic solutions, whose properties are 
discussed in [611Z11I3], arise on setting oq = 0. 

In m we have proven, for any value of the adS radius of curvature i, the existence 
of dyonic soliton and black hole solutions of the field equations fl2.10af[T10dp in a 
neighbourhood of the trivial (Schwarzschild-)adS space-time fl2.17p . Providing the non¬ 
trivial solution is sufficiently close to the trivial solution, the magnetic gauge field 
function Ci;o(r) will have no zeros. In figures [T] and [2] we show two typical nodeless 
solutions: a soliton and a black hole solution respectively. 
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Figure 1. Typical dyonic soliton solution with £ = 1, ai = 0.08715 and W 2 = 0.2. We 
plot the electric gauge field function ao{r) (blue, solid), magnetic gauge field function 
wo(r) (red, dotted) and metric function So{r) (purple, dashed). Both the electric gauge 
field function ao{r) and the magnetic gauge field function wo(r) are monotonically 
increasing and neither has any zeros for r > 0. The metric function mo(r) is not 
shown - it too is monotonically increasing. 



Figure 2. Typical dyonic black hole solution with £= 1, rh = 1, a'f^ = 0.09974 
and ujh — 0.9. We plot the electric gauge field function ao(r) (blue, solid), magnetic 
gauge field function a;o(r) (red, dotted) and metric function So(r) (purple, dashed). 
The electric gauge field function ao{r) is monotonically increasing and the magnetic 
gauge field function uJo{r) is monotonically decreasing. Neither gauge function has any 
zeros for r > Vh- The function <5'o(r) varies only a little: its value on the horizon is 
0.9974, while at infinity So{r) —>■ 1. The metric function mo{r) is not shown - it is 
monotonically increasing as r increases. 
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More detailed properties of the space of dyonic solutions of su(2) EYM in adS can be 
found in [6],|Tlj. Our focus in this paper is the nodeless dyonic solitons and black holes. In 
the case of purely magnetic solutions, it has been proven [HIE] that at least some nodeless 
solitons and black holes are stable under linear, spherically symmetric, perturbations 
of the metric and gauge held. We therefore expect that at least some nodeless dyonic 
solutions will also be stable under linear, spherically symmetric, perturbations. In the 
next section we derive the equations governing such perturbations before proving the 
existence of stable dyonic solutions in section 01 


3. Perturbation equations 

We now derive the equations satished by linear, time-dependent, spherically symmetric 
perturbations of the static equilibrium dyonic solitons and black holes discussed in the 
previous section. 


3.1. Linearized perturbation equations 


We begin with the time-dependent held equations fl2.2af[2Y7j) for the metric (12.5p and 
gauge potential fl2.8p . The appropriateness of the form fl2.5p for studying perturbations 
of the static equilibrium solutions is discussed in section 13.21 The Einstein equations 
(12.2 gp are 

..2 . .2 1 


m = 


252 


a 


'' (7 + a) ] + ( 7 ' + jdy 


m = 2 /i \uuj’ + uy + a) ( 7 ' - 1 - /5)] , 


^ _ 

~S ~ r \: 


00 '^+ 00^ a+ /3) 


+ 


rp2^2 


[u/ -t- (7 + af] , 


and the Yang-Mills equations fl2.2^p take the form 

0 = (^a" - /3'j -h ^ 2 rp - ^ - /3^ - (a + 7 ), 


0 = r^ (fi — a] — 


r^S 


(3-a')+2fiSW (/3 + 7 '), 

{fiSy 


0 = - OJ {y + a) - 2u {y + a) + -—+ a) + {fiS) u ( 7 " -F (3') 

fib 

+ [2 {fiSf u' + fiS {fisy u] ( 7 ' + y ), 

0 = — a) -I- J Cj + uj{a + yy + uj” + fiS {fiS)' uj' 

fib 

- f^S'^oj {y + 7 ')^ (1 - ca. 


(3.1a) 

iy.ib) 

(3.1c) 

(3. Id) 
(3.1e) 

(3.1/) 

(3.1^) 


where a dot' denotes partial diherentiation with respect to time t and a prime ' denotes 
partial diherentiation with respect to the radial co-ordinate r. 
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The ansatz for the gauge potential A^. fl2.8p possesses a residual su(2) Lie algebra 
gauge freedom. If g(t, r) is a diagonal 2x2 matrix depending on t and r, then the 
following gauge transformation leaves the form of the YM gauge potential (12.Sh invariant, 
but changes the matrices A, B and C\ 

+ B ^ B + Q 

C-C^ ^ 0-^ {C - C^) 0, C + C^ ^ 0-^ {C + C^) 0. (3.2) 

Under this gauge transformation the YM gauge held strength Frv transforms as 

Fru Q~^F^^q. (3.3) 

In studies of the stability of purely magnetic su(iV) EYM solitons and black holes 
(see, for example, [ID]), the matrix A is identically equal to zero for the equilibrium 
solutions, and in this case the residual gauge freedom fl3.2|) is used to set M = 0 for the 
time-dependent perturbations as well. Such a choice of gauge simplihes the analysis of 
the resulting perturbation equations in that case. 

Here we are interested in dyonic equilibrium solutions for which the matrix A 
does not vanish. In this case the most appropriate choice of Lie algebra gauge is 
not immediately apparent. Instead of choosing a gauge, we consider gauge-invariant 
variables which do not change under the Lie algebra gauge transformation (13.21) . The 
gauge-invariant variables we use are: 

= ^ = 7 +a, ?7 = 7'-f/3, (3.4) 

together with the variable u which is unchanged by the gauge transformation (13.2p . We 
can also eliminate the gauge-invariant variable since 

'4’ = e-V- (3.5) 

With these new variables, the Einstein equations (l3.1all3Tcll take the more compact 
form 

+ ( 3 - 6 «) 

m = 2yL (uu'+ , {3.6b) 

^ = + + (3.6e) 

and the Yang-Mills equations (l3.1df[3T^ also simplify: 

0 = A ii' -v)- (3.6<0 

W S , 2yiS^uj^ 

0 = - V + ^ - ^ ^ -V) -^— V, 3.6e 

b r‘^ 

0 = — ujf — 2uf + (/i*S')^ {(jjrj' + 2u'ri) + uf + yiS {fiS)' ut], (3.6/) 

yib 
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0 = — a) H- —bj + + (/iS')^ bj” + yiS iyS)' bj' — (yS)'^ bbrj^ 

yb 

+ ■ (3-6^) 

We now consider linearized perturbations about the static equilibrium solutions 
discussed in section [221 The held variables are written as sums of the time-independent 
equilibrium quantities (denoted by a subscript 0, for example yoi'^)) and small time- 
dependent perturbations (denoted by a 6, for example Sy{t,r)). For the static 
equilibrium solutions, we have /So = 7 o = 0 and hence, for the gauge-invariant variables 
(El]), we have 

fo{r) = ao(r), rio{r) = 0. (3.7) 

The time-dependent held variables are therefore written as follows 


/i(f, r) = /io(r) -h 5/i(t, r), S{t, r) = S'o(r) -F 5S{t, r), 

f{t, r) = ao(r) + 6f{t, r), 7]{t, r) = 6r]{t, r), 

bj{t,r) = bJo{r) + 6u{t,r). (3.8) 

The held variables (I3.8|l are substituted into the held equations (l3.6af[3T^, working 


to hrst order in the perturbations. The resulting equations are simplihed using the 
equilibrium held equations fl2.10afEl0dj) - The linearized perturbed Einstein equations 
are then 


6y’ = 


1 f2a 


2/ ,2 

O'^O 


V 




2u'f - 1 

AbJo 


'q -1] 6y+ ^3 


r/io 


+ ra'^ 1 


1-^0 


OCn 




i fo- 


2raQ 

C2 

*^0 


(5^' - 6fj) 


rpo'S'o 

4po 




6fi = - [uq6u -F aou^dr]) , 

AanU? 


5S'= - 

+ 


•■0^0 
rylSo 
4aon;o 
r/igS'o 


2, ,2 






l2 






ryf^Si 


0 


and the perturbed Yang-Mills equations take the form 


0 = {se - s,)') + - 


r^yo 


h0«0ec' 1 /^Otto'S'o ‘iO'oUo 

ob + 


+ MO (Y - Sfi) 


So 


SI 


’-6S- 


0= -6ii + 6i' - ^6S - ,57], 

ho r2 


0 = bOodf + 2Q;o<5d; — aoUo 


6y ^ 6S 
ho ho 


- y^SQUodr]' 


(3.9a) 

{3.9b) 

(3.9c) 


6u 


(3.9d) 

(3.9e) 
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0 


/io5'o \2^oSqu:'q + (/io5'o)^Ci;o] 5t], 

6u! + ^Sq6uj" + fioSo (/io>S'o) duj' + ^qSqujq i^SoSfi' + ^q6S') 


s^o M + 


Wo (1 - 1 :^ 0 ) , 2a§wo 

t^oSl 


+ 


5yL 


f 2aoa;o 

V -s-o 


+ iiqSqUJq 


+ On + 


^iqSI (1 - 


5oj + 2aoWo5^. 


(3.9/) 

5S 

(3.9^) 


3.2. Space-time diffeomorphism gauge transformations 


In addition to gauge invariance with respect to representations of the su(2) Lie algebra 
as discussed above, we must also consider the question of gauge invariance with respect 
to space-time diffeomorphisms. This is relevant when dealing with perturbations of 
a space-time, as we must ensure that the quantities encountered are indeed genuine 
perturbations, and not just artefacts of an inhnitesimal co-ordinate transformation 
carried out on the background space-time. The key idea is to consider the effect of 
such an inhnitesimal co-ordinate transformation generated by the vector held such 
that 


x" ^ x" + 1/L 


(3.10) 


Under such a transformation, the metric perturbation Spru changes according to 


Sg TU Sg TV + ^v9 rW) (3.11) 

where is the Lie derivative along V. Likewise, the gauge potential perturbation (an 
su(2)-valued one-form) undergoes the transformation 


6Aj- — 6Aj- -|- EyA.^. 


(3.12) 


We are considering time-dependent, spherically symmetric perturbations of the 
static, equilibrium conhgurations. Unlike the situation that holds for non-spherical 
perturbations (decomposed into multipoles of appropriate valence), we do not have a 
complete set of gauge-invariant quantities to work with [15]. We must therefore be 
cautious in identifying genuine perturbations, and perturbations which are pure gauge. 
We must also ensure that we are considering perturbations of maximal generality. In 
the co-ordinates (t,r), the most general bare spherically symmetric perturbation of the 
metric has the form 


dgru 


(3.13) 


( Sgoo Sgoi Q \ 

^goi ^gii 

Q Sg22 0 
\ ^ 0 55(22 sin^^ 

where O 2 is the 2x2 zero matrix. Under the gauge transformation (I3.1ip generated by 




(3.14) 


we have 


^goi —t 5501 — EoSqx' + Po ^y, 
dg 22 -t dg 22 + 2ry. 


(3.15a) 

(3.156) 










Stability of dyons in EYM 


11 


We exploit these to simplify the form of the perturbation as follows. We begin with 
the completely general perturbation fl3.13p . We then make a diffeomorphism gauge 
transformation generated by (13.141) . choosing 
^922 


y = 


2r 


(3.16) 


This brings us to a gauge in which 


6922 = 0 . 

This condition is preserved by further gauge transformations provided that y 
fl3.14p . We then apply such a further transformation, choosing V so that 

/ _ 

“ yoSl 

This yields a gauge in which 


(3.17) 
= 0 in 


(3.18) 


^9m = 5<J22 ^ 0, (3.19) 

which is preserved by further gauge transformations generated by gauge vectors of the 
form 


W = (a;(f), 0,0,0). (3.20) 

This represents the generator of the only gauge freedom that remains in the problem, 
and corresponds to a redehnition of the time co-ordinate via fl3.10p . We will refer to the 
gauge condition fl3.19p as the diagonal gauge, and we note that the perturbed metric 
now has the form 


5gru = Diag((55(oo, dgn, 0, 0). (3.21) 

Thus the metric perturbation may be represented by a perturbation of the background 
metric functions: 


/io(r)/ro(r) -h (5/r(t,r), So{r) ^ S^ir) + 5S{t,r), (3.22) 

as assumed in the previous subsection. 

For perturbations of static, spherically symmetric space-times, the choice of 
diagonal gauge (and the reasons behind this choice) are standard and well-known. 
However, for our analysis it is important to understand the residual diffeomorphism 
gauge freedom given by fl3.20p . in particular its effect on the metric and matter 
perturbations. These are summarized in the following lemma. 

Lemma 1. The most general spherically symmetric perturbation of the metric 
may be written in the diagonal gauge h3.21\) where 

^9oo = — 2poSo6S, dgn = ‘^dp. (3.23) 

Under the remaining gauge freedom of infinitesimal co-ordinate transformations 
generated by U3.2(A) . the metric perturbation functions transform as 

dp —)■ dp, dS ^ dS Sqx. (3.24a) 

Furthermore, we can deduce from U3.1‘A] that the matter perturbations transform as 

doo —)■ doj, df) ^ df) -\- aQX, df ^ df + aox, drj —)■ dp. (3.246) 
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The results of this lemma, in particular the behaviour of 5S fl3.24ap and fl3.24&p . 
will be useful in the next subsection for understanding the hnal form of the linearized 
perturbation equations. 


3.3. Linearized perturbation equations in standard form 


We now seek to set the linearized perturbation equations fl3.9af3.9g' ) into a form 
amenable to proving the existence of stable equilibrium solutions. From fl3.9M 13.9 d 3.9j^ 
it can be seen that the perturbation drj is out of phase with the other perturbations: 
when this quantity appears in the perturbation equations with an even (respectively, 
odd) number of time derivatives, all other variables appear with an odd (respectively, 
even) number of time derivatives. (In physical terms, this means that if dr] follows a 
sine wave, the other variables follow a cosine wave, and vice versa.) We therefore dehne 
a new quantity 6k hy 

6k(t, r) 


6r]{t,r) = 


(3.25) 


S-oSo^o 

We note that 6k is dehned only up to an arbitrary function of the radial co-ordinate 
r, and this freedom in dehning 6k will be useful in our later analysis. The equilibrium 
functions of r are introduced in (13.251) because they will enable us to ultimately set the 
perturbation equations into a standard form. With the substitution fl3.25p . equations 
fl3.9f)l I3.9el 3.9yp then take the form 

4 

6ja = — —pp (^ij.qSquJq6uj -f 0!qujq 6ki ), 


0 = 


rSo 
6 k 


fXoSoUo r 




0 


= 6k' + ^6k + ^ 
^0 L'oSq 



2ar] 


■6u — 


Uo 


L'qSo ho'S'o 




(3.26a) 

{3.26b) 

(3.26c) 


Integrating fl3.26ap with respect to time gives the metric perturbation 6jj.{t,r) to be 

4 


6fi = 


rSr 


{fiQSoUQdu -f OiQLijQSK) -f 6J~{r), 


(3.27) 


where 66F{r) is an arbitrary function of r. For purely magnetic background solutions 
with tto = 0, the expression fl3.27l) reduces to that in [10] for this metric perturbation. 

Substituting for (5/i(f, r) (13.271) in (13. 26 cl) and integrating with respect to time gives 
the other metric perturbation 6S{t,r): 


aoUo V aou;^ ) 


+ 2S'o (-f 

Wo r 


1 


Sl 

aQ 


6f -f 6Q{r), 


(3.28) 


where 6Q{r) is another arbitrary function of r. We note that this expression is not valid 
for purely magnetic background solutions with ao = 0. When the background solutions 
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are purely magnetic, the perturbed Einstein equations can be used to find an expression 
for 5S' but not 5S [TO] . 

Substituting for SS{t,r) in (13.266|1 and integrating with respect to time gives the 
following: 


0 = hfj 


{fipSofa'^ 

ap 


6k' — jopSpuipSf' + ‘IfipSpOiQ I 1 + 


2upUj'q \ 

r J 


Soj 


tioSpa'pUp ^^ ^ ( 2[ipSlJl ^ Aapa'^ujl _ ^ 

ap \ r'^ r ao^o J 

(3.29) 


where 6'H{r) is a third arbitrary function of r. A second equation involving 6k can be 
derived from fl3.9ap . substituting in for 6f] using fl3.25p . for h/i using fl3.27p and for 6S 
using fl3.28p . Subtracting the resulting equation from fl3.29p gives a constraint on the 
arbitrary functions 6iF{r), 6Q{r) and 6'H{r): 


0 = 6J^' 


2ao^o 


1 - 2uj, 


/2 


r 


2ra'Q 

T^p 


hoQ^o + 


2ao^o 




6g 


+ 


2ra'n 


-6H. 


A second, independent, constraint on the functions 66F{r), 6Q{r) and 6'H{r) is 
derived from fl3.9(ip as follows. First rearrange the equation resulting from integrating 
fl3.26&P with respect to time to give an expression for 6S, which involves both 6Q and 
61-1. Next differentiate this with respect to r and then substitute into fl3.9dp . simplifying 
using the forms (13.27113.28P of the metric perturbations. Using fl3.29p to eliminate 6k 
from the resulting equation gives the constraint, whose most compact form reads 

1 


^ 2SpapUQ 

^ 9 9 / 

^ /^o«o 

+ 


6g' + ^6g + — ^-^6n' 


Sp 


S'oOiQUp 


fXpapUJp 


k'O 


2S'p 


^ 1 6n. 
^0, 


. (3.31) 

k'O >J0 ^0 / 

Next we use the remaining perturbed Einstein equation (13.9 cp to give a first order 
equation for 6f, which involves only derivatives with respect to r: 


a'n 


0 = 6e-^6f- 

ap 


+ 


^pSp f a[ 


Up 

Aapu'^ [up{l- ul) 


r 

., 2 , 


r'^/j.pu'p 


Up 

1 

+ - + 
r 




S'. 


Up \ ap p-o 

2uq Up Sp 


_o 

^0 


6k' 


+ 


4:apUp 
rjopSp 
rulS^a'ou'o 


4a;2 


+ 


t^oS^p (1 


Un 


4:ralu^ 




+ 


rfJ^pSpUd 

2alu^p 


_ 1 + ^ 

r ap 

., 3 , ,2 


h-Q 

ho 




6K + ^-^6:F + ^6g' 
rfi^oSp Sp 


2ap 


«o^o 


rSp \ulSl 


- w, 


l2 


6g. 


(3.32) 
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We use 03.321) to eliminate from 03.29P and the remaining perturbed Yang-Mills 
equation p.Qg'P . In terms of the usual “tortoise” co-ordinate r*, dehned by 
dr* 1 
dr 1^0 Sq’’ 

the resulting pair of coupled perturbation equations takes the form 
— Sk = 


/\ O y • O 

dl 5 k + 2aodr,5u + Si5k -|- £25uj H- 


rulSo 


-5u) = 


t^oO^o^oOG + - - - - Wo J 

d'^^Su — 2aQdr^5K + S^Sk + S^Su 


5G + 5U, 


fioSWo^J^' + 


2 alojQ 

rSn 


V + 2ci;oWo) 5G 


(3.33) 


(3.34a) 


+ 


^2^0 (1 


w, 


O' + (r + 2n;ow') + p'5^0;' 5T, (3.346) 


r- rpo 

where a prime ' denotes differentiation with respect to the radial co-ordinate r, and the 
coefficients 8 i, 82 , 83 and 84 are given by 


81 — ai + 


lioSl (1 + ul) 2 ^ilSlu'f Aalul ( 2 a' 


+ 


+ 






ao 




82 — 2 po*S'o«o 


ASoaou'^ (1 - u^) 


So. 

(3.35a) 


83 

8 , 


AnoSoaoUoUj'Q _ 1 _ r _^ 

r V“o r 2 ujI /xq So) ’ 

82 — 2 /io*S'oao, 

2 h'O'S'o (1 — Swo) dpgS'oWo^ /po S'q\ 

3“o-^-^ (- + ^ j 


(3.356) 

(3.35c) 


- [2a;o (l - Wg) + r/igWg] . (3.35d) 

For purely magnetic background solutions with ao = 0, we have £^2 = 0 = £^3 and the 
equations fl3.34al [3.346^ are no longer coupled, instead giving separate equations for 5k 
and 5u. In this case 81 reduces to the potential in [7] for sphaleronic sector perturbations 
of the purely magnetic background solutions, while 8 ^ reduces to the potential in [7] 
for gravitational sector perturbations of the purely magnetic background solutions. For 
black hole solutions, the boundary conditions fl2.13p ensure that £^j, z = 1,.. . 4 all vanish 
as r —)■ r/i. For soliton solutions, as r —)■ 0 the functions 82 and 83 tend to constants; 
however 81 and £”4 diverge like a positive constant multiplied by r“^. At infinity, they 
all tend to constants. The functions £^j, z = 1,.. .4 are all regular for r > 0 or r > r/i, 
as applicable, provided that the equilibrium magnetic gauge held function Wo(r) has no 
zeros. Our stability analysis in the next section will be applicable only to equilibrium 
solutions for which this is the case. 

Returning to the case of dyonic background solutions with nontrivial ao, we have 
two coupled, dynamical, perturbation equations (I3.34a[ I3.346P for the perturbations 5k 
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and duj, together with a constraint eqnation fl3.32p for which does not contain any 
time derivatives. We also have in onr system three arbitrary fnnctions of the radial 
co-ordinate r only, namely Q and Ti, which are constrained by the eqnations (I3.3n[ 
I3.3ip . We now argue that all three of these can, without loss of generality, be set equal 
to zero. 

We start by noting that the variable 6k is dehned by fl3.25p only up to an arbitrary 
function of r. This freedom in the dehnition of 6k can be used, via fl3.27p . to set 6 iF = 0. 
The hrst constraint equation (I3.30p then gives 6 % in terms of 6 Q: 


6H = 


‘2cyniijQ . 


r^a 


0 


6g. 


(3.36) 


Substituting for <577 from fl3.36p into the second constraint equation fl3.3ip then gives a 
hrst order differential equation for 6 Q: 


o = 6g'+{--^- 


3S!, 3af, 2uj' 


P-o So 

which can be readily integrated to give 

got^oSle^ 


+ 




+ 


Uo 




r^/ioCTo 


6g{r) = 


'■^^0^0 


where is an arbitrary constant and 

2ao(r')a;o(r')^ 


X = 


'=ro r'2/io(r')a(,(r') 


dr'. 


(3.37) 


(3.38) 


(3.39) 


The lower limit tq in fl3.39p . and the consequences for 6 g, depend on whether we are 
considering soliton or black hole solutions. 

For black hole solutions, the expansions (I2.13p mean that the integrand in (I3.39p is 
regular as r —>■ In this case we set the lower limit of integration tq to be and then 
X = 0(r — Th) as r ^ Vh- At the black hole event horizon, we require 6 fi = 0 and 6 u to 
be hnite: we note by Lemma 1 that this statement is gauge invariant. Therefore, from 
(I3.27p . recalling that we have already set 6 iF = 0, it must be the case that {r — rh)6K —)■ 0 
as r ^ r/i, so that both 6 k and (r — rh)6K' are integrable at the horizon. We also 
require the perturbations 6 S, 6'if and 6 f to be hnite at the event horizon. While these 
quantities change under an inhnitesimal co-ordinate transformation generated by (I3.20p . 
from fiTM [TMI we can see that such a transformation maintains the hniteness of 
these perturbations at the horizon. Considering the hrst YM perturbation equation 
fl3.9dp as an ODE in r for 6 ijj = 6f' — 6ri, carrying out a formal integration, and imposing 
the condition that 6 'ijj remains hnite at the horizon, we hnd that (r — rh)~^6f must be 
integrable at r = Vh- Therefore, from (I3.28p . we see that 6 g must also be integrable at 
the horizon. However, from the dehnition fl3.38p . the properties of X as r —)■ r/j and the 
boundary conditions (I2.13p . we see that 

2 goSi 


6g = 


Th) r^a'^ 


+ 0 ( 1 ) 


(3.40) 


as r —)• r^. Thus 6 g fails to be integrable at the horizon unless = 0- Therefore both 
6 g and 577 must vanish identically in the black hole case. 
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For soliton solutions, from the expansions fl2.12p we see that the integrand in 03.391) 
is 0{r~^) as r —)■ 0. In this case we choose the lower limit of integration to be tq = 1. 
With this choice, as r —)■ 0 we have X = 2 logr + (9(l) and therefore = O(r^) as r 0. 
Substituting this and the expansions 02.121) into 03.381) . it can be seen that 6Q = 0{r~^) 
as r —)■ 0 unless = 0. In order to keep the origin regular, the perturbation 5uj must 
remain hnite at r = 0 , and we must also have 5/r—)- 0 asr—)- 0 to avoid a curvature 
singularity. As in the black hole case, these are diffeomorphism-invariant statements. 
From 03.27P and the boundary conditions 02.121) we see that dn also remains hnite as 

r —)■ 0, and, as a consequence, rdn' —)■ 0 as r ^ 0. We also require 55', dfj and to 

be hnite at the origin. As in the black hole case, using the results of Lemma 1, this 
requirement does not change if an inhnitesimal co-ordinate transformation generated by 
03.2np is applied. With this assumption and the above behaviour of 5 k and its derivative 
as r —)■ 0, equation 03.281) then implies that, at worst, r‘^5Q —?• 0 as r ^ 0. Therefore, 
as in the black hole case, the only possibility is = 0 , so that both 5Q and dH vanish 
identically. 

Setting ddF = 0, 5^ = 0 and dV. = 0, the perturbation equations 03.34al 13.3451) 
simplify to 

— dk = — d‘^^dK + 2aodr,duj + Sidn + S2duj, (3.41a) 

—did = — d'^^du — 2aodr,dK + S^dn + Sidu, (3.415) 


which do not involve the perturbation df. Once perturbations 5 k, du solving 03.41 of 
13.4151) have been found, the perturbation df is computed by solving the constraint 
equation 03.321) . which now simplihes to 
1 


tto 


-5^ 


= (5^(5k, 5a;) 
hO'S'o 


-5k" 


-5a;' 


aoUo 


Uo 


ho “So 


f U p (1 - a;g) ^ 1 ^ 
r^fipUjQ r 


0^0 

ap 

r 


ho 

ho 


or 

’^0 




2uq lip Sp 


4aoa;o f r fipS^ (1 - o;^) 


rfipSp 
rfilS^a'ou'o 

Aalul 


4a;2 


+ 


ho 

1 


Araluil 


r ao 


5k' 


duj 


ho 


Sk 


ho 


0 

5-0 


+ 


2alut 


5k, 


(3.42) 


where we have dehned a new quantity (^(5 k, 5a;) which depends on 5 k, du and the 
equilibrium solutions. 

The fact that we do not have a dynamical equation for df can be understood from 
the analysis of section 13.21 As discussed in that section, we have a residual inhnitesimal 
diffeomorphism gauge freedom generated by O3.20|) . which corresponds to a redehnition 
of the time co-ordinate. Such a gauge transformation changes df according to fl3.245p . 
Since this is a gauge transformation, the perturbation equations should be independent 
of the choice of the gauge function x{t). When df changes as in fl3.245p . the quantity 
(ag ^5.^)^ does not change, and so the equation (13.421) remains invariant. If we had a 
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dynamical equation for involving time derivatives of such an equation would not 
be invariant under the gauge transformation fl3.246j) . Since we have already eliminated 
Sf), the only other perturbation which changes under the residual diffeomorphism gauge 
transformation is dS fl3.24ap . However, the changes in (13.24 and 5S fl3.24ap mean 
that the equation fl3.28p for 5S remains invariant under the gauge transformation. 

Our purpose in this paper is to prove the existence of dyonic equilibrium soliton 
and black hole solutions of the static EYM held equations which are stable under linear 
perturbations satisfying (13.41 al 13.41 I3.42p . We turn to this proof in the next section. 


V = 


8 = 



(4.2) 


4. Existence of stable dyonic solitons and black holes 
f.l. General argument 

The perturbation equations (13.41 al 13.41 M can be written in the compact form 

— V =—dl^v+ Vdr,v + 8v =Uv, (4.1) 

where v = {6 k, doj)^ is the vector of perturbations, and the matrices T) and 8 are given 
by 

0 2Q!q 

—2oo 0 

with 8 i,... ,84 given in fl3.35af[T35dp . We restrict our attention to static equilibrium 
solutions for which the magnetic gauge held function Ci;o(r) has no zeros. In this case 
the functions 81 ,... ,84 are regular for all r > 0 in the soliton case and all r > r^ in the 
black hole case. They diverge as r —)■ 0 for soliton equilibrium solutions, but this is not 
an issue, as discussed after fl4.39p below. 

In (14. ip we have dehned an operator W by 

G = -d^^+Vdr^+ 8 . (4.3) 

The operator Id will be symmetric if 

= -V and 8 ^ = 8 - dr,V. (4.4) 

The hrst of these conditions is clearly satished. With the form of P in (14.2p . the second 
condition is satished if 

82 — 8-4 = dr, (2q!o) = 2/ro‘S'oaQ, (4.5) 

which can be seen to be the case from (I3.35cp . Therefore W is a symmetric operator. 

To derive sufficient conditions for 7/ to be a positive operator, we hrst write it in 
the form 

u = x\-sv, 

where the operator V does not contain any derivatives, 

X = Sr. - Z, (4.7) 


(4.6) 
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Z = 



Zl2 

0 


We choose the functions Z 12 and 2^21 such that 


(4.8) 


dr,Zi2 — S2, dr^Z2l — 8^. 


(4.9) 


Then the matrix V takes the form 

0 

0 8 i-Z 



(4.10) 


Therefore the operator lA will be a positive symmetric operator if the functions /Ci and 
/C 2 are positive everywhere, where 

/Cl = - ^2^, /C 2 = ^4 - ^12- (4.11) 


If W is a positive symmetric operator then the static equilibrium solutions will be stable 
under linear, spherically symmetric perturbations. It turns out that there are three 
related stability properties: these will be described in detail in section 14.41 

For the moment, let us assume that to prove the stability (as will be characterized in 
section 14.41) of the static equilibrium solutions it is sufficient to prove that the functions 
/Cl, K ,2 (14.111) are positive everywhere and, in addition, have the asymptotic and other 
properties used in the stability arguments in section 031 Before proving the existence of 
dyonic solutions for which this is the case, we present a couple of numerical examples, 
considering the static equilibrium solutions shown in hgures [1] and [2l Having found an 
equilibrium solution by integrating the static held equations (l2.10af[T10dl) . the hrst step 
in showing that they are stable is to numerically integrate (14.9p to hnd the functions Z 12 
and Z 21 . In hgures [3] and 0] we show the results of these integrations for the equilibrium 
solutions presented in hgures [T] and |2] respectively. Note that the equations (14. 9 p only 
dehne the functions ^12 and Z 21 up to the addition of an arbitrary constant. This 
constant is chosen so that the functions vanish at either the origin, or the black hole 
event horizon, as applicable. 

Once we have the functions Z 12 and Z 21 , we then compute the quantities /Ci and /C 2 
(14.lip . If we can show that /Ci and /C 2 are positive for all r, then the operator U (14. 6 p is 
a positive symmetric operator and hence, as argued above, the equilibrium solutions are 
stable. In hgures [5] and [6] respectively we present the results of calculating /Ci and /C 2 
for the dyonic soliton and black hole solutions shown in hgures [1] and |2l In both hgures 
Eland El it can be seen that the two functions /Ci and /C 2 are positive everywhere. 

We now turn to the proof of the existence of dyonic soliton and black hole solutions 
which are stable under linear, spherically symmetric, perturbations of the metric and 
gauge held. The stability criterion that we apply is that the operator lA governing 
the evolution of the perturbations (14.61) is a symmetric, positive operator. As we have 
seen, this condition reduces to establishing that the functions /Ci and /C 2 are positive 
everywhere - that is, for all r > in the black hole case, and for all r > 0 in the soliton 
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Figure 3. Functions Z 12 (blue, solid) and Z 21 (red, dotted) defined by (14.91) for a 
typical nodeless dyonic soliton solution with .^ = 1, ai = 0.08715 and UJ 2 = 0.2. The 
additive constant of integration has been chosen so that both functions vanish at the 
origin. For this example, Z 12 is monotonically increasing with r; on the other hand 
2^21 has a maximum and is negative for sufficiently large r. 



Figure 4. Functions Z 12 (blue, solid) and Z 21 (red, dotted) defined by (14.9|) for 
a typical nodeless dyonic black hole solution with £ = 1, r/j = 1, = 0.09974 

and ujh = 0.9. The additive constant of integration has been chosen so that both 
functions vanish on the horizon. As in the example shown in figure [U in this case Z 12 
is monotonically increasing as r increases. The other function, -2^21: has a minimum 
just outside the event horizon and then slowly increases as r increases. 
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Figure 5. Functions /Ci (blue, solid) and IC 2 (red, dotted) defined by (14.111) for a 
typical nodeless dyonic soliton solution with £ = 1, oi = 0.08715 and W 2 = 0.2. Both 
functions are positive for all r. For small r, they exhibit very similar behaviour, only 
becoming distinguishable for r ~ 1. As r —^ 0, both functions diverge, in accordance 
with (|4.35|) . They both have a minimum, and then tend to (positive) constants as 
r —>■ 00 . 



Figure 6. Functions /Ci (blue, solid) and IC 2 (red, dotted) defined by (14.111) for a 
typical nodeless dyonic black hole solution with £ = 1, m = I, a'j^ = 0.09974 and 
ujh = 0.9. Both functions are positive for all r > r^- Both vanish on the horizon, then 
rapidly increase to maximum values, before tending to (positive) constants as r —>■ 00 . 
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case. We establish this in each of these cases by recalling some of the results from m, 
in which we proved the existence of the background space-times being studied here. We 
consider the black hole case hrst. 


4-2. Existence of stable dyonic black holes 


In m, we proved the existence of static, spherically symmetric, asymptotically adS 
black hole solutions of the su(2) EYM equations for which the metric and gauge held 
functions /io. So, ao and cjq have the following properties: 

(i) {fio,So,ao,uJo) e C^{[rh,+oo),'Mf) x +oo),M^); 

(ii) /io(’"/i) = 0 and ^o{r) > 0, a;o(?") > 0 and So{r) > 0 for all r > Vh', 

(hi) there is a continuous mapping of the initial data (imposed at the horizon) 
{Sh,uJh,o:'ff) G to the solution (/xq, (So, ao, i^o) £ C'^([’"/i,+oo), M^) x 

C^{[rh,+oo),R'^)-, 

(iv) this continuous mapping yields a unique global solution of the held equations in a 
neighbourhood of the trivial datum {Sh,uJh,o:'h) = (1)1)0) G which gives rise 
to Schwarzschild-adS space-time as an embedded solution of the system. We will 
refer to this as the trivial solution: it nevertheless plays an important role in what 
follows. 


The trivial solution is characterized by fl2.17p . so that 


/ro(r) = 1 - 


2M A , 

-r 

r 3 


A simple calculation using fl3.35af[T35djl then yields £2 = £z = ^ and 


£i — £i — — [1 — 


2M A . 

- E 

r 3 

By a choice of constant of integration, we can then take Z 12 = Z 21 = 0, and so 


/Cl = /C2 = -7 1 


2M A . 

- E 

r 3 


(4.12) 


(4.13) 


(4.14) 


It is then immediate that both /Ci and /C 2 are positive outside the horizon. 

For the non-trivial background solutions, we can use the asymptotic forms 02.131) 
to establish that 


Si(r) = 0(r - rk), r-> r*. i = l,2,3,4, 


(4.15) 


and 


£i{r) = £i^oo + 0{r r ^+ 00 , i = l,2,3,4. 


Then 


and so we can take 


= 0(1), r^Th, / = 2,3, 


£2{r' 


£3{r') 


f,„(r')S„(r')*'’ X.„,Ao(r')S„(r')*'’ 


(4.16) 

(4.17) 


(4.18) 















Stability of dyons in EYM 


22 


giving 


Zi 2 {r) = 0{r-Th), Z 2 i{r) = 0{r-Vh), r ^ th, (4.19) 


and 


Zi 2 {r) = 0{r ^), Z 2 i{r) = 0{r ^), r-)■+cx). (4.20) 

It follows that, as r ^ Th, 


ICi{r) = Si{r)+ 0{[r-Th]^), K: 2 {r) = S^^r) + 0{[r - Vhf), (4.21) 

and, at infinity. 


fCi{r) = £i{r) + 0{r ^), /C 2 (r) = £^ 4 (r)+0(r ^), r +co.(4.22) 


Therefore /Cj, i = 1, 2 are continuous as r ^ r^,. It follows from fl3.35af[T35dl I4.18P and 


from points (i) and (ii) in the list above that the /C* are continuous on (r/j, +cxo). Thus 
ICieC^[rh,+oo), i = l,2. (4.23) 


(4.23) 


Furthermore, it follows from this observation and from point (iii) above that for each 
Th > 0, the /Cj depend continuously on the initial data {Sh,uJh,o:'if) G Thus we can 
conclude the following: 

Proposition 2. For each A < 0 and Vh > 0, there exists a neighbourhood U of the trivial 
initial data point {Sh,uJh,o:'ff) = (1,1,0) G such that for all {Sh,uJh,o:'ff) G U, the 
corresponding unigue black hole solution of the static EYM eguations, whose existence 
is guaranteed by Proposition 7 of ms. satisfies 


/Cj(r) > 0 for all r > Vh, * = 1, 2, 


(4.24) 


and thus is linearly stable under spherically symmetric perturbations. 

From (I4.15( I4.19p . the functions /C*, i = 1,2 vanish at the event horizon. We note 
that a necessary condition for them to be positive outside the horizon is that 



(4.25) 


Using the expressions fl3.35al 13. 35 dp . we find that 



+rl{l - Arl){3ujl - 1 )] } • 


(4.26) 


Expanding fl4.26p about the trivial data shows how a neighbourhood of allowed values 
of the initial data for which ki, i = 1,2 are positive may arise: 



(4.27) 
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Positivity of the /Cj in the soliton case follows by a similar argument. 

In [H], we proved the existence of static, spherically symmetric, asymptotically 
adS, soliton solutions of the su(2) EYM equations for which the metric and gauge held 
functions po, Sq, ao and uq have the following properties: 


(i) (po) >S'o, ao, Wo) G +oo), M^) x Ycxo), M^); 


(ii) fio{r) > 0, a;o(r) > 0 and S'o(r) > 0 for all r > 0 and these functions satisfy the 
asymptotic behaviour of fl2.12l) at the origin; 

(iii) there is a continuous mapping of the initial data (imposed at the origin) 

(S'i,W 2 ,ai) G to the solution (po, *S'o, ao, wo) G C^([0, +cxd), M^) x 

C'2([0,+cx)),M2); 

(iv) this continuous mapping yields a unique global solution of the held equations in a 
neighbourhood of the trivial datum (S'i,W 2 ,ai) = (1,0,0) G which gives rise to 
adS space-time as an embedded solution of the system. This is the trivial solution 
in the solitonic case. 

For the trivial solution, we hnd £^2 = ^^3 = 0 and 



(4.28) 


Then as in the black hole case, we can take = ^21 = 0 in the trivial background, 
and positivity of /Ci = and /C 2 = on (0, -|-cxd) is immediate. 

For the non-trivial background space-times, we can use (12.121 I3.35M I3.35cp to 
establish that 


E 2 {r) = 2aiSi{l + 4 :U 2 ) + 0{r), E^i^r) = SaiSiu ;2 + 0{r), (4.29) 


as r —)■ 0. The general form of the functions E 2 fl3.356p . E 3 f|3.35cp and the properties 


of the background metric functions listed above, along with the limiting behaviour at 
the origin, then shows that E 2 and E 3 are continuous on [0, -|-cxd). The same is true of 
EijifioSo) for i = 2,3 (again appealing to (I2.12p and the properties listed above) and so, 
from (14.9p we can take 



(4.30) 


which yields 


- 2q:i(1 + 4 ^ 2 )?’+ o(r^), r0. 


(4.31) 


Likewise, 



(4.32) 


with 


^21 (^") = 8 aiW 2 r-I-O(r^), r —)■ 0. 


(4.33) 
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We note also that the asymptotic behaviour at inhnity as described by fl2.16p yields 

Zu = 0{r-^), Z 2 i = 0(r-i), r ^+cx). (4.34) 

At the origin, and Z 21 are dominated by Si, £^ 4 , and we hnd that 
2 92 2 9^ 

JCi = + 0(1), A ^2 = + 0(1), r —)■ 0, (4.35) 

and so the fCi, i = 1 , 2 are clearly positive in a (one-sided) neighbourhood of the origin. 
At inhnity, we hnd 

/Cl = 0(1), /C 2 = 0(l), r^+ 00 . (4.36) 

With these details in place, positivity of the /Cj, z = 1, 2 on (0, -|-oo) follows in this 
case in exactly the same way as the black hole case, giving: 

Proposition 3. For each A < 0, there exists a neighbourhood U of the trivial initial data 
point (S'i,a; 2 ,«i) = (1,0,0) G such that for all (S'i,a; 2 ,ai) G U, the corresponding 
unigue soliton solution of the static EYM eguations, whose existence is guaranteed by 
Proposition sofun, satisfies 

/Cj(r) > 0 for all r > 0, i = 1 , 2 , (4.37) 

and thus is linearly stable under spherically symmetric perturbations. 

4.4- Nature of the stability 

In the previous section we have proven the existence of dyonic soliton and black hole 
solutions of su(2) EYM such that the functions /Ci and /C 2 fl4.1ip are positive everywhere. 
Therefore the operator U fl4.6l) is a positive symmetric operator. In this case we have 
three related stability properties, which we now discuss in detail. 

4.4- 1. Mode stability. First, we consider time-periodic perturbations with frequency a, 
so that 

5K{t,r) = 6uj{t,r) = e*'^*(5a;(r), (4.38) 

(we return to Sf shortly). Then the perturbation equations fl4.ip for v = {6n{r), 6u{r))'^ 
take the form of a standard Schrodinger-like eigenvalue problem: 

Uv = a'^v. (4.39) 

To establish that all the eigenvalues must be positive, we need to impose suitable 
boundary conditions on the perturbations. We set 6K{r) and 6u{r) to vanish at the 
origin (for soliton solutions), event horizon (for black hole solutions) and at inhnity. 
With these conditions on Sn^r) and Su{r), the inner products {v,v) and {v,Uv) are 
hnite. As usual, the inner products are dehned by integrals over the range of the 
tortoise co-ordinate r*: for solitons this is r* G [0,r*max] for some constant r*max while 
for black holes it is r* G (—oo,0]. The only issue that may be of concern is the form 
of the zero-order terms in (14. ip in the soliton case, which diverge as r~^ at the origin 
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fl4.35p . However the coefficient of is sufficiently large and positive to ensure that 
solutions vanishing at the origin decay rapidly enough to yield integrable terms in the 
inner products. 

With suitable boundary conditions in place, the positivity of the operator U implies 
that all the eigenvalues are positive, so that a is real and the perturbations do 
not grow exponentially with time. There is one further remaining subtlety, namely 
that we do not have a dynamical equation for 6 f{t,r). Suppose however that we have 
perturbations Suit, r), Su(t, r) of the form (I4.38p which satisfy (14.391) with > 0. Then, 
by integrating fl3.42p . we find 

pr 

= oio{r)X{t) + ao{r) / ^{6K,6u)dr', (4.40) 

J r'=ro 

where X{t) is an arbitrary function of time t and the quantity ^{6 k, 6u) depends linearly 
on the perturbations {6k, 6u) and can be found in full in fl3.42p . Since the perturbations 
{Sk,6oj) are time-periodic (I4.38p . and given the form of ^{6k,6u;) in (I3.42p . the second 
term in fl4.40p does not grow exponentially with time. The hrst term corresponds to 
an inhnitesimal diffeomorphism generated by (I3.2np . and by a suitable choice of x{t) in 
fl3.24&D can be transformed away. Therefore 6 f{t,r) is also time-periodic with > 0 
and does not grow exponentially with time. The fact that the metric perturbations 6 yi 
and SS also do not grow exponentially with time can be deduced from their explicit 
forms (13.27113.28P in terms of the matter perturbations. 


4.4-3- Finite energy. Next, we note that the form (14.61) of the perturbation equations, 
with V positive, gives rise to a positive, conserved energy for the system. That is, 
dehning 

E[v]{t) = ^{v,v)+ ^{v,Uv), (4.41) 

we have 

E[v]{t) > 0 for all v, 

= 0 r; = 0. (4.42) 

Positivity follows from the decomposition (14.6p . the dehnition of the adjoint operator 
yl and positivity of V: 

E[v]{t) = ^||h||2 + ^Wxvf + ^{v, Vv), (4.43) 

where the norm is that associated with the inner product. We note that the factor | 
appearing in (I4.4ip is conventional, to make the link with the conserved energy of a 
particle in classical mechanics. Conservation of the energy corresponds to the fact that 

4 {i5M(()} = 0. (4.44) 

The derivation of this equation is straightforward, relying on the symmetry of V in 
(H.lOp . and on the identity 

{v, X^Xv) = ixv, X'i’) = ^ {{xv, xv) + {xv, xv)) ■ 


(4.45) 
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E[v\{t) = i?[i>](0) =: Eq for all t>0. 


(4.46) 


Thus a second stability property holds: the system gives rise to a positive definite energy, 
which is conserved by the evolution. In particular, no blow-up of the energy is possible. 

4.4-3. Pointwise bounds. Conservation of the positive dehnite energy can be used to 
show that the perturbation v is bounded throughout the evolution. The argument that 
follows applies to those background solutions described by Propositions 2 and 3 above. 
These solutions give rise to functions Z 12 , ^ 21,^1 and /C 2 with the properties required 
for the arguments below. 


It follows from fl4.43p and fl4.46l) that 

\\dr^v — Zv\\‘^ < 2Eo, {v,Vv)<2Eo. 


(4.47) 


We note that Eq is the energy of the initial conhguration and that we have used (14.7p . We 


can deduce a pointwise stability result using these bounds. Slightly different arguments 
are required in the soliton and black hole cases. 

Solitons. In the case of solitonic dyons, we show that (I4.47p enables us to obtain an 


a priori bound on the El^ norm of v. Then Sobolev embedding leads to a pointwise 


bound. To see this, we note that 



= Ci(n, Vn) 
< 2CiEq, 


(4.48) 


where 


Cl = max{ sup /C^ ^(r), i = 1, 2}, 


(4.49) 


re[0,oo) 

whose existence and positivity is guaranteed by the argument of Section 4.3 above, where 
we show that the fCi are positive and continuous for r G (0, -|-cxd), and satisfy (14.351 IT36ll . 
Note that the /Cj depend only on the background solution functions. Similarly, 


Zv\\ < C2||n||, 


(4.50) 


where 


C 2 = max{ sup |Za(r)|, a = 12, 21}. 


(4.51) 


rG[0,oo) 
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Again, the argument of Section 4.3 shows that this term is hnite - see 04.301 - l4.34p . 
Then 

= \\dr,v - Zv + Zv\\ 

< \\dr^V — Zv\\ + \\Zv\\ 

< \/2Eq + C2\/^^/‘2Eq, (4.52) 

where we have used 04.47114.481147501) . Thus we have an a priori bound for the El^ norm 
||n||i of v: 


||n||i := \\dr,vf + ||n||^ < 2C^Eo, 

(4.63) 

where the positive constant C 3 , determined by Ci and C 2 , depends only on the 
background solution. Then the Sobolev inequality |T 6 ] 


(4.54) 

yields a pointwise bound for the perturbation v: 


\v{t,r)\ < C 3 \/^ for alH > 0 ,r G [ 0 ,+cxo). 

(4.55) 

We note that in these inequalities, - is the Euclidean norm in 
expresses a pointwise bound on the soliton perturbation v. 

and so fl4.55p 

Black holes. A slightly diherent argument is required in the black hole case, due to the 
diherent behaviour of the functions /Cj, i = 1, 2 at the horizon. The general mathematical 
argument is the same, relying on the Sobolev inequality to derive a pointwise bound 
from a bound on the norm. First, we note that 

1 < \\9r.,'V — Zv\\ + \\Zv\\ 

< ^/Wo + C^{{v,Vv)f/^ 

< (1 + C 4 ) \/ 2i7o, 

(4.56) 

where 


C 4 = max < sup , sup — > , 

I r£[rh,oo) TG[r’^,oo) ^2 J 

(4.57) 


which is positive and hnite by virtue of the properties listed in 04.15l - IT.22p . Of particular 
importance is positivity of the /Cj, and the hnite limits of the relevant ratios at both the 
horizon r = r/j and at inhnity. 

Next, we introduce the positive matrix 

^/^ 0 
0 VKf 


W = = 


(4.58) 


and we dehne 


w = Wv. 


(4.59) 
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= {v,Vv) < 2 Eq. 


(4.60) 


Also, 


where 


and 


< C4dr,v\\ +<^ 6 ^ 11 , 


C 5 = max sup \/^, sup 

r£[rh,oo) rG[r^,oo) 


Cq = max < sup 

( re[r^,oo) 

1 

= - max sup 

rG[rh,oo) 






, sup 

re[rh,oo) 


dr,'/}^ 




hO'S'o 


/c;(r) 


/Cl 


, sup 

r£[rh,oo) 


hO'S'o 


/C^(r) 


/Co 


(4.61) 


(4.62) 


(4.63) 


both of which are positive and hnite by virtue of (14.151 - 14.221) . Then (14.56114.60114.6ip 
establish an a priori bound for the norm of w, and the Sobolev inequality (14.541) 
applied to w yields 

\Wv\<Cj^/¥o, (4.64) 


where the constant C 7 is constructed from C 4 — Ce, and depends only on the background 
solution functions. Let us recall certain properties of the /Cj established in Section 4.2 
above: these functions are positive and continuous on r G {rh,+oo), with hnite limits 
as r —)■ +CXD. They vanish at the horizon, with the asymptotic behaviour tCi = 0{r — rh) 
as r ^ rh- It follows from (I4.58P and (I4.64p that v is bounded on any interval of the 
form [rh + e, Tcxo) with e > 0. The boundary condition lim^^r^ v = 0 ensures that this 
pointwise bound extends to the entire interval [r/j,+oo). This establishes a pointwise 
bound for v{t, r) for all f > 0 . 


5. Conclusions 

In this paper we have proven the existence of dyonic soliton and black hole solutions 
of four-dimensional su(2) Einstein-Yang-Mills theory in asymptotically anti-de Sitter 
space which are stable under linear, spherically symmetric, perturbations of the metric 
and gauge held. Although the static, dyonic, equilibrium solutions of the held equations 
were found numerically over hfteen years ago [ 6 ], their stability has not been investigated 
until now. 

The perturbation equations for linear, spherically symmetric, perturbations of the 
dyonic equilibrium solutions are much more complicated than those for purely magnetic 
equilibrium solutions, which is perhaps why the stability question has not been explored 
previously in the literature. In the purely magnetic case, with a suitable choice of 
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Lie algebra gauge, the perturbation equations decouple into two sectors, known as 
the “sphaleronic” (odd-parity) and “gravitational” (even-parity) sectors |T^. This 
decoupling of the equations for odd- and even-parity perturbations greatly simplihes 
the analysis (see, for example, [5] for the su(2) case and [TUI for the larger 5u{N) 
gauge group). In the purely magnetic case this decoupling arises because the static 
equilibrium solutions are both spherically symmetric and invariant under a parity 
reflection x —)■ —x of the space co-ordinates and the purely magnetic gauge held. 
Therefore the perturbations which are odd and even under the parity rehection can 
be considered separately. While the dyonic static solutions considered here are still 
spherically symmetric, the gauge held is not invariant under a parity transformation 
(see, for example, the discussion in section 2 of [3]) and so the perturbation equations 
for the odd- and even-parity perturbations no longer decouple. 

We have therefore taken an alternative approach in this paper, and considered 
perturbations which are invariant under Lie algebra gauge transformations. The 
perturbations of the metric functions can be found in terms of the gauge held 
perturbations and hence eliminated. After some manipulation, the equations for the 
remaining three gauge held perturbations can be cast into a pair of coupled Schrodinger- 
like equations fl4.39p involving just two of the perturbations {Sn and 6u) and a single 
constraint equation fl3.42p which does not contain any time derivatives and determines 
the third perturbation (5,^) once the other two are known. The lack of a dynamical 
equation for is due to a remaining diheomorphism gauge freedom, corresponding to 
a redehnition of the time co-ordinate. 

In m we proved the existence of static, spherically symmetric, dyonic soliton and 
black hole solutions of the su(2) EYM equations in adS for which the single magnetic 
gauge held function uq has no zeros. These nodeless solutions exist for any value of the 
negative cosmological constant A < 0, in a neighbourhood of the trivial (Schwarzschild- 
adS) embedded solution. By analysing the afore-mentioned pair of coupled Schrodinger- 
like perturbation equations, we have been able to prove that nodeless dyonic solutions 
sufficiently close to the embedded trivial solution are stable under linear, spherically 
symmetric, perturbations of the metric and gauge held. This extends the proof of the 
existence of stable purely magnetic solutions to the dyonic case. 

It would be interesting the explore the consequences of the existence of stable dyonic 
black holes for the “no-hair” conjecture, in the form stated by Bizon [18], namely 

Within a given matter model, a stable stationary black hole is uniquely 

determined by global charges. 

Combining the results of [ini|19], we have evidence that the above conjecture is true for 
purely magnetic black holes in 5u{N) EYM in adS. To investigate whether the stable 
dyonic black holes, whose existence we have proven here, satisfy the above conjecture, 
appropriate electric and magnetic charges would need to be dehned, and then one would 
need to determine whether these charges uniquely characterize the black holes. We hope 
that this question will be the subject of further research. 
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Dyonic black hole solutions of EYM in adS have received a great deal of attention 
recently in the literature when the event horizon, rather than being topologically 
spherical as considered in this paper, has planar topology. In [20] it was found that 
there is a second order phase transition from the embedded planar Reissner-Nordstrom- 
adS solution to a black hole with a non-trivial dyonic YM condensate. There is now 
a substantial literature on such planar dyonic EYM black holes as models of p-wave 
holographic superconductors (see, for example, |2l] for a selection of references, and 
[22] for a recent review). The thermodynamic behaviour of these planar dyonic EYM 
black holes has been studied extensively (this is the key aspect of their interpretation 
as holographic superconductors) but little is known about their classical stability. At 
least some purely magnetic topological black holes in su(2) [23] and sxi{N) [2l] EYM 
in adS have been proven to be stable when the magnetic gauge field functions have no 
zeros. In this paper we have considered the stability of spherically symmetric dyonic 
solutions of 5u(2) EYM, and a natural question would be to extend this to topological 
black holes, or to a larger gauge group (the existence of dyonic solitons and black holes 
in su(Y) EYM in adS having recently been established [25]). However, our results in 
this paper are for dyonic solutions where the magnetic gauge held function ujq has no 
zeros, whereas for the solutions of relevance for holographic superconductors ujq has a 
single zero, located on the adS boundary. This will complicate the classical stability 
analysis and therefore we leave this as an open question for future investigation. 
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